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ANALYSIS OF SHAPE OF POROUS COOLED MEDIUM FOR AN 
IMPOSED SURFACE HEAT FLUX AND TEMPERATURE 
by Robert Siegel 
Lewis Research Center 

SUMMARY 

A two-dimensional porous cooled medium has part of its boundary open to a coolant 
reservoir and another part subjected to uniform heating by an external source. Coolant 
from the reservoir is forced through the medium and exits through the part of the bound- 
ary subjected to the heat load. For design purposes it is specified that the coolant exit 
surface is to be at a given temperature. This analysis provides a method to obtain the 
shape of the porous surface to meet these simultaneous constraints of a specified tem - 
perature and heat loading. The analysis includes temperature variations of fluid density, 
fluid viscosity, and matrix thermal conductivity. The solution is obtained by combining 
the energy equation with Darcy's law in such a way that a potential function can be intro- 
duced that satisfies Laplace's equation. All of the heat -transfer and coolant flow quan- 
tities can be expressed in terms of this potential, which is governed by simple boundary 
conditions. The determination of the shape of the porous cooled region is thereby re- 
duced to a free -boundary problem such as in irrotational, inviscid free jet theory. Two 
illustrative examples are carried out: a porous leading edge with coolant supplied 
through a slot and a porous cooled duct with a rectangular outer boundary . 


INTRODUCTION 

Some of the advanced power producing devices such as fusion reactors and gaseous 
nuclear reactors will require advanced heat -transfer designs to cope with the large wall 
heat fluxes that are anticipated. One proposed cooling scheme is transpiration or porous 
wall cooling. A metal or ceramic wall is made in a porous form and coolant is forced 
through the wall so that it exits from the side of the wall exposed to the high heat flux. 
Any heat conducted into the wall is transferred to the coolant and carried back out. The 
result is an effective cooling technique. If part of the wall heat flux is provided by heat 



transfer from a convective boundary layer, the transplant blowing out of the wall tends 
to push the boundary layer away from the wall, thus reducing the convective load. 

The power required to pump the coolant through the wall represents a loss. Also, if 
the hot gas heating the porous wall is being used as the working fluid in a power produc- 
ing device, the mixing of the coolant with the gas decreases the gas temperature and thus 
reduces the thermal efficiency of the device. In a device that is not ground based, such 
as a high speed aircraft, the coolant must be carried along, and this represents an addi- 
tional weight penalty. These considerations illustrate that the coolant should be used as 
effectively as possible. A poor distribution of coolant flow will tend to overcool some 
portions of the surface and undercool others. 

This report will consider a two-dimensional porous cooled region. The objective of 
the analysis will be to determine the shape of the region so that it will provide a specified 
heat -transfer performance. Since the region is two-dimensional, its length in the direc- 
tion normal to the cross section must be large in comparison with a characteristic di- 
mension of the cross section. A porous region along the leading edge of a wing or turbine 
blade would be an example of this type of configuration. The coolant enters the porous 
region through the part of its boundary that is open to a reservoir. The coolant exit sur- 
face has a heat flux specified along it. It is desired to maintain the coolant exit surface 
at a specified temperature, for example, at a value slightly below the melting point or 
one dictated by structural strength considerations. A method is presented here to obtain 
the shape of the porous surface such that a given surface temperature is maintained while 
the surface is being heated uniformly by an external source. 

The analysis includes the effect of variable properties. The fluid viscosity can be a 
function of temperature, and, if the coolant is a gas, its density is found from the perfect 

gas law. The thermal conductivity of the solid matrix material can also be a function of 
temperature. 

To obtain an analytical solution for a situation including all of these effects, some 
simplifying assumptions are also required. The coolant mass flow through the porous 
material is assumed to be at a pore Reynolds number small enough so that Darcy’s law 
applies. The thermal resistance between the fluid and the porous matrix material is as- 
sumed small enough so that the local fluid and matrix temperatures are equal. As a re- 
sult of this assumption a single energy equation can be written for the balance of energy 
transport by conduction in the matrix and by convection of the flowing coolant. 

The variable fluid properties cause the flow and energy equations to be coupled. 
Because the coolant exit boundary of the porous medium has an unknown shape, a solution 
must be found that will simultaneously satisfy the flow and energy equations and at the 
same time determine the unknown boundary shape so that it is compatible with the spec- 
ified boundary conditions. The solution is obtained by use of a suitably defined potential 
function that satisfies Laplace's equation. The temperature distribution and local heat 
flux can be directly related to this potential. This was shown in reference 1 for a porous 
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region of specified shape. The analysis of reference 1 is extended here to include var- 
iable thermal conductivity of the matrix material and to account for the unknown shape of 
the coolant exit surface. 

The shape of the unknown boundary is obtained by using conformal mapping such as 
in inviscid free jet theory. The author previously adapted the free-jet theory to deter- 
mine the free -boundary configuration of a two-dimensional freezing layer (ref. 2). It 
will be shown that the results of the freezing problem can be directly carried over to the 
porous medium problem. To illustrate the application of the method, results are carried 
out for a two-dimensional leading edge with coolant introduced through part of its bound- 
ary and for a porous cooled duct with a rectangular outer boundary. 


ANALYSIS 

Geometry and Specified Conditions at Boundary 

The general type of porous region being considered is shown in figure 1. The con- 
figuration is two-dimensional, that is, it is a cylinder of arbitrary cross section, infi- 
nitely long in the z direction, and with its axis parallel to that direction. A coolant res- 
ervoir at pressure and temperature p^jt^ is adjacent to a portion Sq of the porous 
material that is open to the flow. Coolant flows through the porous medium and leaves 
through the surface S, which is exposed to a heat flux q g . The remaining surface area 
is impervious to the flow and is also insulated. Some specific geometries of this type 
will be discussed later (see figs. 3(a) and 7(a)) to illustrate the application of the 
analysis. 

The surface S through which the coolant exits is specified as being at uniform tem- 
perature and also as having a uniform heat flux being transferred to it. This would ap- 
pear to be an overspecified set of boundary conditions, but it is the shape of S that is 
not specified. The objective of the analysis is to obtain the shape of S that will allow 
the simultaneous specification of these two conditions. This would provide the optimum 
shape to accommodate the incident heat flux for the specified temperature. For example, 
a heat flux could be specified and the shape found that would cause the surface temper- 
ature to be within a given value below the melting point or be at a temperature fixed by 
structural considerations. If the porous material can sublime without disturbing the 
porosity, the shape would be that formed by sublimation as a result of being subjected to 
the specified heat loading. In this instance the t g would be the sublimation temperature. 

As the fluid in the reservoir at p,^ approaches the coolant inlet surface Sq, the ac- 
celeration effects are neglected with regard to affecting the fluid pressure. Hence, the 
pressure along S Q is equated to the reservoir value, p Q = p^. Since p^ > p g , the cool- 
ant flows from Sq to S. (Symbols are defined in the appendix. ) Inside the porous 
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medium, the fluid picks up the energy being conducted in at the surface S and then car- 
ries it back out of the medium. Since p Q and p g are both constants, the fluid veloc- 
ities along both the inlet and exit surfaces are locally perpendicular to these surfaces. 
This fact enters the boundary conditions along these surfaces. 

The effective thermal conductivity of the matrix material is k which is a function 
of temperature and is based on the entire cross-sectional area rather than on only the 
area of the solid material. Because the conductivity of the matrix material is generally 
much larger than that of the coolant, the heat conduction in the coolant is neglected. The 
Darcy velocity u, that will be used throughout the analysis is the local volume flow di- 
vided by the entire cross-sectional area rather than by the pore cross-sectional area. 


Governing Equations 

The size of the porous medium is usually quite large compared with the size of the 
individual pores, so that it is not necessary to deal with the individual pores and volume 
averaged continuum equations can be used. Then the following equations apply within the 
porous material for compressible and incompressible coolants: 

Conservation of mass - 


v • (pu) = o (compressible) 
V • u = 0 (incompressible) 
Darcy's law (fluid inertia effects are neglected) - 

u = - Vp 
M(t) 


(la) 

(lb) 


( 2 ) 


Conservation of energy (kinetic energy is neglected) - 

V • q = 0 (3) 


where 


q = -k m (t) Vt + puCpt (4) 

Perfect gas law (for compressible case) - 

P = pR t (5) 
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Boundary Conditions 


As the fluid in the reservoir approaches the porous medium, it receives any heat 
that is being conducted from the coolant inlet face Sq back into the fluid. This raises 
the coolant temperature (and hence matrix temperature, since they are locally equal in 
this analysis) to tg, which is an unknown quantity along the surface Sq and which will 
be obtained later in the analysis and will be found to be independent of position along Sq. 
The thermal conductivity of the fluid is generally much smaller than that of the porous 
matrix material, and hence, as the fluid approaches the wall, the temperature rise from 
t^ to tg takes place in a thin layer compared with the thickness of the porous material. 
A boundary condition is obtained by balancing the heat conducted out of the surface with 
the energy carried back into the matrix by the fluid. The pressure along the boundary 
can be taken at the reservoir value since the pressure change from flow acceleration in 
the reservoir is small compared with the pressure drop through the porous material. 
Hence, the conditions along the coolant inlet face are 


k m S 0 ’ * = P C p (t “ *->0 ' U 


P = P 0 = Poo = constant 


for (x,y) on S Q 


( 6 ) 


On the surface S of unknown shape and through which the fluid leaves the porous 
material, the temperature and imposed heat flux are both specified as constant so that 
the boundary conditions are 


t = t g = constant ^ 
k m rig • Vt = q g = constant > 
p = p = constant 

On the impervious surfaces the conditions are 



for (x,y) on S 


( 7 ) 


for (x,y) on 


(8) 



Equations and Boundary Conditions in Dimensionless Form 


The governing equations and boundary conditions will now be placed in dimensionless 
form. Note that the functions P and M differ for the compressible and incompressible 
flow cases; these differences will result in the two cases reducing to the same set of 
relations. The dimensionless variables are 



Y 


y_ 

h_ 


T 



P 


incompressible 


compressible 


(9) 


M(T) = < 


r \ Lt(t) ^°° . ... 

— — incompressible 

2 Uoo t 

p(t) 

Moo 

V. 


compressible 


(10) 


k (t) 
K (T ) = JH11 

m' 


( 11 ) 


qh 

Q = r 


k„ t 


(12a) 
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(12b) 


Q = 


q b 


s r 


k t 
m. oo oo 


PooCn^oo 
\ = Z 

^Poo^m, oo 

v/ii + i^ 

ax aY 


(13) 


For the compressible case equation (5) is first used to eliminate p from equa- 
tions (la) and (4). Then for both the compressible and incompressible cases, Darcy's 
law is used to eliminate u. The result is that equations (1), (3), and (4) reduce to 


= _L VMT • VP 

(14) 

MT 


V • Q = 0 

(15) 

(Vi r^J 

■K VT - ~ VP 
m ■»» 

(16) 


M 


By a similar manipulation boundary conditions (6) to (8) become 


K T A ~ \ a ~ 

-2— n n • VT + — n n • VP = 0 
T - 1 U M u 


> for (X, Y) on S. 


P = 1 


J 


0 


T=Tg A 




K m n s • ^ = Qs 


> for (X,Y) on S 


P = P 


s J 


l, • VP = O' 


n t • Q = 0 


for (X,Y) on S. 


(17) 


(18) 


(19) 
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Formulation in Terms of Potential Function 


At this point we shall postulate, pending later verification that P is a function of T 
only. It will be found that with this condition a solution is obtained satisfying the govern- 
ing equations and boundary conditions, and hence this is the required solution. It will 
now be shown that this assumption permits the energy flux to be expressed as the neg- 
ative gradient of a potential $ , 



(20) 


and then from equation ( 15) 


V 


2 


$ 


= 0 


( 21 ) 


Hence the potential function is a solution to Laplace's equation; the boundary conditions 
for <1> will be given later . 

Since P is only a function of T 


1 

M(T) 


VP = 


dP ~ 
— VT 


M(T) dT 


Hence, if we put 


it follows that 


/ T 
u 

/*T 

g(T)= / 


1 dP 
M(T) dT 


dT 


K m (T)dT 


Vh(T) = VT = — L_ ^ VT = L_ VP 

dT M(T) dT M(T) 


Vg(T) = VT = K (T)VT 
dT m 


(22a) 


(22b) 


(23a) 


(23b) 
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Equation (16) can then be written as 


/sj 


Q = -Vg(T) -XVh(T) = -V[g(T) + Xh(T)] = -V<f- 


(24) 


where the potential function $ is defined by 


4> - 4? 0 = g(T) + Xh(T) = 


n 


K m (T) + — ~ 
m M(T) dT 


dT 


(25) 


and <f>Q is an arbitrary constant that can be used to fix the level of the potential. Since 
P is only a function of T, is a function of T only, and hence P and T are func- 
tions of 4> only. 

Equations (15) and (16) have been expressed in terms of a potential as given by equa- 
tions (21) and (20) with the potential given by equation (25). Now equation (14) will be 
considered. This can be written as 


v 2 p - 1 d < MT) v* • dp v* = |v*i 2 

MT d$ d$ MT d4> d$ 




Now V Z P = V • VP = V • — V<f> . In view of equation (21) this becomes" 
’ d$ 


V 2 P=^|v*| 2 

d* 2 


(26) 


(27) 


From equation (25) 


dT m M dT m M d$ dT 


(28) 


Equations (27) and (28) are used to eliminate P from equation (26) giving the following 
equation, which is a restatement of equation (14) in a form that determines T as a func- 
tion of 4>, 


d 2 T 1 dT 

| 

L/k - T d M dT l 

j.2 K T d4> 
a<t> m 

\ m dT / d$ 


*Note that V 


=^v. v$+v$. v® = ^ v 2 <t> V$ 

d$ d$ d$ d$ d ^,2 


^v 2 * + ^]v$| 2 


d$ 


d$ 
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or equivalently 


— I (k 

d<f> T \ m d$ 



( 29 ) 


Equation (29) can be integrated to obtain T($), but before doing this the first bound- 
ary condition of equation (17) will be expressed in terms of It follows from equa- 
tions (24) and (23) that 


— VP = V$ - K VT 
M m 


Substituting into the boundary condition gives 


K m T . ~ . ~ 

— n 0 ' ST+V (V4-K m VT)=0 

/V( 

Using the relation VT = (dT/d4>) V<i> results in a new expression of the boundary 
condition: 


K 


1 + 


m dT\ - ~ 

'n Q -V^=0 for (X,Y) on S Q 


T - 1 d4>; 


(30) 


Solution for Temperature Ratio in Terms of Potential Function 

Equation (29) can be integrated twice to obtain 


4 > 



K m dT 
1 + CjT 


(31) 


The constants C ^ and are to be evaluated from the boundary conditions. Thus in 
serting the solution (31) into equation (30) gives 


(1 + C l ) n Q - = 0 for (X, Y) on S Q 
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This boundary condition is satisfied if C ^ = -1 so that equation (31) gives the relation 
between T and $ 


4 > 



( 32 ) 


Relation Between Pressure Function and Temperature Ratio 


The second of boundary conditions (17) involves P. It is necessary to relate P to 
T so that P can be related to $ and the boundary condition formulated in terms of $ . 
From equation (28) 


dP = /d$ _ K \ M 
dT VdT m X 


(33) 


and from equation (32) 


d$ _ K m 
dT 1 - T 


(34) 


These relations are combined to obtain 


dP 

dT 


' 1 A K m M , T K m M 
1 - T /X 1 - T X 


(35) 


Integrating this from Sq to an arbitrary position in the medium and taking into account 
the second boundary condition ( 17) gives the relation between P and T , 


P - 



(36) 


The value Tq is unknown in this integral. To evaluate it, the first and third boundary 
conditions of equation (18) are imposed to give 
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T 


(37) 


P 


s 



1 - T 


dT 


This integral can be carried out once the thermal conductivity and viscosity variations 

are specified so that K m (T) and M(T) are known. The quantity Tq is thereby related 

to the specified quantities P and T . 

s s 


Boundary Conditions for Determination of Potential Function 

The temperature has been obtained in terms of a potential $ (eq. (32)), and the 
pressure has been found in terms of the temperature (eqs. (36) and (37)). Hence all 
quantities will be known if the potential can be found as a function of position in the po- 
rous material. To obtain the simplest boundary conditions for the potential <£, which is 
governed by Laplace’s equation, the constant $q in equation (25) can be set equal to 
zero to give 


$ 


i K <T) + 


X dP 
M(T) dT 


dT 


and since T = Tq along Sq the boundary condition is 

$ = 0 for (X,Y) on Sq 
W ith 4 >q = 0 equation (32) becomes 


<f> = 


X i-i 


dT 


This is evaluated at the coolant exit face to obtain the boundary condition 


* = * s = 


f T s K m 

X 1 -T 


dT for (X, Y) on S 


(38) 


(39) 


(40) 


(41) 
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The surface temperature T is specified and Tq is calculated from equation (37). 
Hence, $ is known. 

It follows from equation (20) and the second condition in equation (19) that $ must 
satisfy the boundary condition on the insulated impervious surfaces 

£. • V4> = 0 for (X, Y) on S i (42) 

Since P is a function only of T, and T is a function only of 4>, the first boundary con- 
dition (19) is then automatically satisfied. 

The remaining boundary condition to be considered is that the free boundary must 
have a shape such that the second condition of equation (18) will be satisfied. The con- 
dition can be written as 


K n • ® V$ = Q 
m s d * s 


(X, Y) on S 


From equation (34) (dT/d$) e = (1 - T)/K so the boundary condition becomes 

s m _ 


v<f> = 


Q c 


1 - T, 


(X, Y) on S 


(43) 


Equations (39) and (41) to (43) provide the required boundary conditions to solve 
equation (21) for <t>. In the solution of Laplace's equation for the potential function, it 
would be convenient to have the boundary potentials on Sq and S go from 0 to 1 rather 
than from 0 to <f> , and to have the normal derivative on the coolant exit boundary be 
unity. To accomplish this, the potential is normalized as 


<p(X, Y) 


4>(X, Y) 


(44) 


Also let Q g /(1 - T s )$ s = D and 


r>u 



Then cp is obtained from 



0 


(45) 
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with the boundary conditions 


<p = 0 for (X, Y) on Sq 

(46a) 

r—i 

II 

9- 


(46b) 

I 

> for (X, Y) on S 


n • V<v? = 1 
s ^ J 

1 

(46c) 

/V 

n. . v (p = 0 

for (X, Y) on S i 

(46d) 


These conditions are summarized in figure 1(b). 

By eliminating T Q from equations (40) and (4 1) the temperature distribution in 
terms of cp is given by 


<P - 1 = 



(47) 


Mass Flux at Coolant Exit Surface 

A quantity of practical interest is the local coolant flux distribution leaving the po- 
rous medium, 


pu " n = - — — pV p 
S M(t) 


n 


for (x, y) on S 


(48) 


First consider the case of a compressible flow. By using equation (5) and introducing 
dimensionless variables this becomes 


pu • £ = - JL. P- Vp • 5 = — nv» _X_ ~p . 5 
p(t) Rt s C p h r MT s 


(49) 


But upon noting that VP = (dP/dT)(dT/d4>)V4> we find from equations (35) and (34) that 


VP = 


KM 
m i 


T \d-T)^ = MT~^ 


vl - T/ K. 


m 
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Inserting this in equation (49) gives 


h C ~ 

r P p u • n = -n ■ V$ for (X, Y) on S (50) 

, s s 

K m , °° 


By use of equation (43) the dimensionless mass flux is given by 


h r C p 


P U c 



(51) 


or 



C p (t s - 


(52) 


A similar manipulation for incompressible flow shows that equation (52) also applies for 
that case. 


Determination of Unknown Coolant Exit Boundary 

The analysis up to this point has shown that the temperature distribution in the po- 
rous medium can be expressed in terms of a potential where the potential is a solution to 
Laplace’s equation. The difficulty in solving the Laplace equation in the porous medium 
is that the coolant exit boundary has an unknown shape. In this problem the heat flux and 
temperature are both given along the coolant exit boundary. The relation between these 
quantities is what is usually desired in a porous medium heat -transfer analysis. The de- 
tailed temperature distribution in the medium is generally of lesser importance but is 
needed in thermal-stress calculations. Hence, the main result to be obtained here is the 
shape of the unknown coolant exit boundary - this type of problem is often called a free- 
boundary problem . 

The solution of free -boundary problems has its foundation in inviscid free-jet theory 
where the velocity potential in the jet flow region is governed by Laplace’s equation. In 
reference 2 a free-boundary theory was developed for use in freezing problems. The un- 
known shape was determined of a steady -state two-dimensional frozen region formed by 
cooling the frozen region along part of its boundary while subjecting another part to a 
convective heating condition. This method was extended to other frozen shapes in 


l 
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reference 3 and to a radiative heating condition in reference 4. (In ref. 4 the specified 
heat flux at the surface is nonuniform. ) 

The present solution has been shown to depend on a potential with boundary condi- 
tions given in figure 1(b). A comparison with the solution method as described in ref- 
erence 2 reveals that these are the identical boundary conditions in the free -boundary 
freezing problem. Hence, this previously developed technique, and the results obtained 
with it, can be directly applied to the present problem. The general analytical technique 
will be briefly reviewed here. Then some examples will be given of porous configura- 
tions that can be applied in practice. 

The analytical technique is based on the use of conformal mapping. Following ref- 
erence 2 the negative potential is taken as the real part of an analytic function: 


W=-(p+i\p (53) 

The lines of constant cp and constant ip form a two-dimensional orthogonal net. By us- 
ing conformal transformations between the W plane (to be described) and the physical 
plane, the transformed functions will always be analytic. Hence in the physical plane the 
potential will satisfy Laplace's equation. Since a function that is analytic in a given re- 
gion is completely determined by its boundary values, it is only necessary to require that 
the mapping satisfy the boundary conditions. 

The potential W plane is shown in figure 2, and the porous region occupies a rec- 
tangle. The coolant inlet and outlet surfaces Sq and S are lines of constant <p 
(eqs. (46a) and (46b)). Along the impermeable S i surfaces the normal derivative of <p 
is zero (eq. (46d)); as a result of this orthogonality, these boundaries are lines of con- 
stant ip. The height of the rectangle is unknown and will be related to the specified heat 
flux along S. If the mapping between the W plane and the physical plane can be obtained 
subject to the constraint in equation (46c), then the line 3-4 in the W plane can be trans- 
formed into the physical plane and the unknown free surface thus determined. 

Let the coordinates in the dimensionless physical plane (fig. 1(b)) be designated by 
the complex variable Z = X + iY. To find the relation between the W and Z planes the 
fact is utilized that the derivative of an analytic function is independent of direction and 
hence the derivative with respect to Z can be written as 


dW _ 3 cp . dip 

dZ 0X 3X 


(54) 


Using the Cauchy -Riemann equation dip/ 3X = 3<p/s¥ and defining a quantity £ gives 


This is integrated to give 


Z = j IdW + C 2 (56) 

where C 2 is an integration constant. If £ can be found as a function of W then the in- 
tegration can be performed to obtain Z which thus relates Z to W. 

To obtain a relation between £ and W a potential derivative plane (-3 cp/dX against 
i d(p/dY) is constructed from the known boundary conditions. In addition to the normal 
derivative boundary conditions shown in figure 1(b), there is the condition that dcp/dl = 0 
along S Q where l is a path along the boundary. This gives along Sq 

d^ = 3^dX + 3^dY = 0 ( 57 ) 

d l 0X dl 0Y dl 

Since Sq has a specified shape, the dX/dZ and dY/d Z are known along this boundary. 

Hence, the curve representing this boundary can be drawn on the £ plane. Along S the 

boundary condition n • Vw = 1 means that this boundary will be part of a unit circle in 

s 

the £ plane. This behavior will be illustrated by a few examples. 

After the representation in the £ plane has been found, the conformal mapping be- 
tween the £ and W planes gives the required £(W) relation so that the integral in equa- 
tion (56) can be evaluated. 

To demonstrate the method and obtain some useful shapes for porous cooled media, 
two configurations will be analyzed: a leading -edge region with fluid supplied through a 
coolant slot and a porous duct with a rectangular outer boundary. 


EXAMPLES ILLUSTRATING APPLICATION OF GENERAL METHOD 
Leading-Edge Region Fed Through Coolant Slot 

This geometry is shown in figure 3(a) and consists of a porous leading-edge region 
with coolant being supplied through a slot in an insulated wall that supports the porous 
region. The coolant exit surface S is to be maintained at uniform temperature t 

s 

under the influence of a uniform imposed heat flux q , which could be supplied, for ex- 
ample, by convection from an external stream. 

The region in dimensionless coordinates is shown in figure 3(b) along with the values 
of cp and its derivatives on the boundaries. From this point on, the analysis exactly 
parallels that for determining the configuration of a frozen layer in reference 2, and 
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these previous results can be used with only a change in nomenclature. The regions in 
the potential and potential derivative planes are given in figure 4 . 

The mapping functions between these planes have been obtained in reference 2, and 
the integration in equation (56) carried out for Z(W). The result was evaluated along the 
coolant exit surface to yield its coordinates as 


:ln 


1 + 0 


Is 2 In 


Vi-o 


2 ,2 


s In 


y/l - 


■ tan 


K b -4 


1 - b‘ 


(58) 


where ^ is a variable with range -b < i; < b. The b depends on the imposed physical 
conditions, and this relation will now be determined. 

In the notation of reference 2 the boundary condition at the free surface was 
3T/3N = 1, and in the present analysis it is n g • V (p = 1. The T and cp are com- 
parable quantities so that the only difference is in the nondimensionalizing factor in N 

ft/ 

compared with that in the V, This shows that the quantity y in reference 2 is to be re- 
placed in the present analysis by h r /D = k m ^(t^ - t g )$ s /q s . Then it directly follows 
from equation (32) in reference 2 that the parameter b is related to the physical quan- 
tities by 


V~g_s -Ql*. 


v 


In 


^1 -b J 


K 


(V ” b 5 ) 


(59) 


Another quantity of interest is the total heat (Qj. 0 j.) being conducted into the surface 
S, which is also equal to the energy transferred to and then carried out of the medium by 
the coolant. This is equal to the imposed heat flux integrated over the area of the free 
surface S. The flow rate per unit area into the medium along S n is given by 

^ j u 

-pu • Hq g . From equation (50) it follows that this can be written for the present geom- 
etry as I 0 


h C pu ■ n n < — * « . 

■ r - P - = n n • V4> =!£ 

V» ° 8Y 


for (X, Y) on S Q 
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To obtain the total flow, integrate over the inlet surface S Q to obtain 


t- 


r c „pu • n Q 
r p r. dx = 2h. 
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(The integration limits refer to the numbers in fig. 3(b). ) By use of the Cauchy -Riemann 
equations this can be integrated to obtain 


/- 


C pu ■ n n 
r P_ °dx = 2$ 


m, °° 


h r°i 

ir i = x 


^ dX = 2$ s h r [<M3) - ^(4)] = 2$ s h r [^(2) - <MD] 


The conformal mapping results (eqs. (24) and (31) of ref. 2) give 


IK2) - Ip(l) = K - i — . 

K ^/jTb 5 ) 


From an overall heat balance 


= - 2 f t P" ' 5 0 C p<*s 


too) d X 


Combining these relations and noting that 4> g is negative, gives 

jtot _ C P g = (60) 

2k m, oo^s ” t °o^4 > s l 2k m,J*sl K(Vl - b 2 ) 

The Qj.^ is per unit length normal to the x-y plane in figure 3(a). The G is the 
total mass flow of coolant per unit time and unit length normal to the x-y plane, that is 

r 3 _ 

G = -2 J pu • n Q dx 

(The negative sign is present because the u is in a direction opposite to n^.) 

For a given set of physical conditions the left side of equation (59) is computed (note 

that the 4> is found from eqs. (41) and (37)). The b can then be found from 
s 
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equation (59), and the coordinates of surface S and the obtained from equa- 

tions (58) and (60). The results are shown in figures 5 and 6. 

Porous Duct with Rectangular Outer Boundary 

As a second example of the conformal mapping technique, results are given for a 
porous cooled rectangular duct (fig. 7(a)). An interior shape of the duct cross section is 
to be obtained that will maintain the interior surface at a given temperature when a uni- 
form heat flux is being transferred to the interior surface by convection or radiation from 
a high -temperature gas. From symmetry, only the lower left quadrant of the duct need 
be considered. 

A comparable free -boundary problem for a solidified region inside a duct was ana- 
lyzed in reference 3. These results can be directly used by only changing the nomen- 
clature as in the previous example. The porous region in dimensionless coordinates is 
shown in figure 7(b). The mappings into the potential and potential derivative planes are 
given in figure 8. For detailed information on these mappings and the conformal trans- 
formation between the £ and W planes the reader is referred to reference 3. A brief 
summary of results is given as follows: 

The aspect ratio of the duct is related to two mapping parameters c and d by 



where K is the complete elliptic integral of the first kind. This is used to find c,d 
pairs that correspond to various aspect ratios. Then the coordinates of the coolant exit 
surface of the porous medium can be found from 




The values of c and d can also be used to obtain the corresponding imposed phys- 
ical conditions from the relation 


k m> cc(t s -tJKJ 


J 2(c + d) g J ( c - l)(d - 1) 

f(c+l)(d+l) L V(c + l)(d+ 1). 

k “(fS) 


(63) 


The total heat transferred to the interior surface of the porous medium (equal to the 
energy carried away by the coolant) per unit length of the duct is given by 
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(64) 


The configurations of the porous region for various aspect ratios are given in fig- 
ure 9. Each of these figures gives curves for various values of the parameter containing 
the imposed physical conditions. Figure 10 gives the total heat transfer and coolant 
mass flow rates as a function of this parameter. 


SUMMARY OF ANALYSIS 

The geometric aspects of the solution are concerned with determining the unknown 
free boundary of a two-dimensional region in which Laplace’s equation is valid (eq. (45)) 
and with boundary conditions (46). This can be accomplished by conformal mapping as 
described in relation to equation (56). 

The mapping requires representing the porous region in the complex potential and 
potential derivative planes. In the complex potential plane the region is always a rec- 
tangle, and the derivative boundary conditions determine its shape in the potential deriv- 
ative plane. The dimensionless heat flow rate as given for the examples by equa- 
tions (60) and (64) involve <& . This is obtained from equation (41) in which the quantity 
Tq is found from equation (37). 

The potential at any location in the porous region can be obtained from the result of 
the integration in equation (56), which relates the coordinates in the potential plane to 
those in the dimensionless physical plane. The temperature at that physical location can 
be found from equation (47), which gives the relation between the temperature and the 
potential. The Tq in this relation is found from equation (37). 
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DISCUSSION 


A method is developed for obtaining the heat -transfer characteristics of a two- 
dimensional porous cooled region having a free boundary. A free -boundary shape is 
found such that it will be at a specified temperature while being subjected to an imposed 
uniform heat flux. The solution was obtained by showing that all the heat -transfer quan- 
tities could be related to a potential function that is a solution to Laplace’s equation in 
the porous region. Conformal mapping could then be used, as developed for free- 
boundary problems, to obtain the configuration of the unknown boundary. The procedure 
can be extended in principle to a nonuniform heat flux at the free boundary. In this in- 
stance the free boundary in the potential derivative plane would not be part of a circle, 
and the mapping into the potential plane would require more sophisticated analytical tech- 
niques or numerical methods. An example of this type of mapping by using a series ex- 
pansion is given in reference 5. In reference 4 a particular nonuniform heat flux is con- 
sidered as imposed by radiant heating. 

Results are given for a porous leading edge region as shown in figure 3(a). This 
type of geometry might be used to cool the stagnation region along the leading edge of an 
airfoil. The shape of the porous region is shown in figure 5 as a function of the param- 
eter k m ^(tg - t 00 )|$ g |/q s a. Figure 5 gives the coolant total mass flow rate and the 
total heat transferred to the coolant exit surface of the porous region per unit length of 
the region. From an overall energy balance the must be transferred to the coolant. 

Hence there is the relation 


Hot = GC p (t s - *»> 

Because of the interaction of the various quantities, the results in figures 5 and 6 
can be discussed in various ways. One approach is to consider a situation where a given 
pressure ratio is available and the surface of the porous medium is not to exceed a cer- 
tain temperature t for a specified heat loading q With these quantities given, the 
parameter k m (t g - t oo )|$ s |/q g a is uniquely determined, and the shape of the porous 
region is found from figure 5. For the purpose of discussion, consider what would hap- 
pen if q g is kept constant and the porous region is made thinner. For simplicity ignore 
property variations so that 4> g remains fixed. Then from figure 5 a thinner region cor- 
responds to a smaller value of the parameter and hence t is decreased. It would ap- 
pear then that an improved cooling situation would result by decreasing the thickness. 
However, it must be realized that according to figure 6 a decrease in the abscissa corre- 
sponds to an increased coolant flow. Hence decreasing the surface temperature to below 
what is really required results in wasted coolant. If the thickness of the porous region is 
made larger than given by the solution, the coolant flow will be decreased and the surface 
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temperature will be higher than desired. Thus the solution gives the porous region shape 
that provides the minimum coolant flow with the correct flow distribution to cope with the 
imposed heat load and not let the surface temperature exceed the temperature limit that 
has been specified for it. 

When the porous region in figure 5 becomes very thick compared with the coolant 
slot width, it approaches a circular shape. This would be expected from the uniform 
boundary conditions. The coolant is then flowing outward in a radial way from what is 
approximately a line source (normal to the x,y plane) at the origin of the region. 

The second configuration treated is a porous cooled duct with rectangular outer 
shape as shown in figure 7(a). The coordinates and boundary conditions in the dimension- 
less system are shown in figure 7(b). The mappings into the potential and potential de- 
rivative planes are given in figure 8. This yields the porous shapes in figure 9 for ducts 
with aspect ratios from 1 to 5. For each aspect ratio a set of curves is given for various 
values of the physical parameter k m ^(tg - t^) s f /q g b . It should be noted that the map- 
ping yields two roots for some values of the parameter. The group corresponding to very 
thick porous regions is shown dashed. The corresponding heat and mass flows are shown 
in figure 10. 

The interesting feature of these results is the fact that for each aspect ratio there is 
a maximum value of the physical parameter beyond which no solution exists. To elab- 
orate on this, consider a situation where t is fixed and q„ is being changed. Fora 

s s 

large q the physical parameter is small, and the solution given by the solid curves in 

b 

figure 9 gives a thin layer. From figure 10 the corresponding flow rate is large. If the 
q g is decreased, the physical parameter is increased. Less coolant is required, and 
the thickness of the porous region is increased to decrease the flow and thus maintain the 
surface at the specified t As q is further decreased, the thickness of the porous 
medium can be increased to reduce the flow and thus maintain the surface at t g . How- 
ever, as the thickness is increased, the inside surface area of the duct is reduced; 
thereby, the total heat loading on the surface is being reduced more rapidly than in direct 
proportion to the reduction of q . A point is reached for a very thick region where to 
maintain a given t the q must be increased with increased region thickness. If q_ 

SS S 

were decreased, the t would decrease below the specified value. This is in the range 

O 

of the dashed curves in figures 9 and 10. 


CONCLUSIONS 

An analytical solution was obtained to determine the optimum shape of a two- 
dimensional porous cooled region. The optimum shape provides the proper distribution 
of coolant flow through the region surface that is subjected to a uniform heat loading in 
order to maintain this surface at a specified uniform temperature. The analysis includes 
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the effects of variable fluid viscosity and density, and it includes the effect of variable 
heat conductivity of the porous matrix material. The solution was found by combining 
the governing equations in such a way that a potential could be defined. This potential 
satisfies Laplace's equation in the porous region and is subject to simple boundary con- 
ditions. All of the heat -transfer quantities can be related to this potential. The potential 
and shape of the porous region are obtained by use of a conformal mapping technique 
analogous to the free streamline method of irrotational inviscid free -jet theory. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 15, 1972, 

502-28. 
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APPENDIX - SYMBOLS 


A,B dimensionless half widths of long and short sides of rectangular duct 

a half width of coolant slot; half width of long side of rectangular duct 

b parameter in mapping; half width of short side of rectangular duct 

Cp specific heat at constant pressure 

C Cg constants of integration 

c,d mapping parameters for rectangular duct 
D dimensionless quantity, Q g /(1 - T g )<i> s 

F elliptic integral of the first kind, F(?j,k) = f drj/^/l - k^ sin^ rj 

* / 0 

G mass flow rate of coolant per unit length normal to x,y plane 

h r reference dimension of porous region 

i,j unit vectors in x and y directions, respectively 

/•ir /2 / / 2 2 

K complete elliptic integral of the first kind, K(k) = I d^/yi - ^ s i n V 

- / 0 

K m thermal conductivity ratio, k m A m 

k effective thermal conductivity of porous region 

m 

l path along boundary S Q in X,Y plane 

M for incompressible case M - (l/2)(jz//i o0 )(t 0O /t); for compressible case 

M = li/V-n 

N dimensionless outward normal coordinate in X,Y plane 


A 

n 

P 

P 

Q s 

S;ot 

Q 

I s 

q 

R 

s 


unit outward normal 

2 

for incompressible case P = p/p*,; for compressible case P = (p/p*,) 
pressure 

dimensionless heat flux imposed at interface q s h r 

heat conducted into porous surface (equal to heat removed by coolant) per unit 
length normal to x,y plane 

dimensionless energy flux vector, qh r A m 
heat flux imposed at porous surface S 
energy flux vector, -k m Vt + puCpt 
perfect gas constant 
coolant exit surface of porous medium 
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8^^ insulated and impervious surface 

S (J coolant inlet surface of porous medium 

T temperature ratio t/t^ 

t absolute temperature 

u velocity vector 

W complex potential -<p + ii// 

X,Y dimensionless coordinates, x/h r , y/h r 

X,Y dimensionless coordinates, XD,YD 

x, y coordinates in physical plane 

Z complex variable X + iY 

£ complex potential derivative, - 4= + i 

SX SY 

k permeability of porous material 

X parameter , Poo C p K Poo /2 Poo k m ^ ^ 

ju fluid viscosity 

k dummy variable in mapping 

p fluid density 

$ potential defined by equation (25) 

$ g potential along coolant exit surface 

cp dimensionless potential <f>/<i> 

s 

x(y function orthogonal to cp 

V dimensionless gradient i _ + j — 

SX SY 

V dimensionless gradient i JL + j JL 

SX SY 

Subscripts: 

i insulated and impervious surface 

s on surface where coolant exits from porous medium 

0 on surface where coolant enters porous medium 

00 at coolant reservoir condition 
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(a) Physical coordinates and imposed conditions. 
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(b) Dimensionless coordinates with boundary conditions in terms of potential. 

Figure L - Two-dimensional porous region with unknown shape of coolant 
exit surface; p^ > p s , t s > t^. 
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Figure 2. - Porous region mapped into potential plane, 

W =-cp+ i ip. 
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(a) Physical coordinates and imposed conditions. 
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(b) Dimensionless coordinates with boundary conditions 
in terms of a potential. 

Figure 3. - Porous leading-edge region with coolant 
supplied through slot. 
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(b) Region in potential derivative plane, £= HdipIdX) + i(dy>/dY). 
Figure 4. - Auxiliary planes used in mapping procedure. 











Figure 6. - Dimensionless heat and coolant mass flow rates at surface of porous 
leading edge region as function of imposed physical parameter. 



(a) Physical coordinates and imposed conditions 



(b) Dimensionless coordinates with boundary conditions in terms of ^ ^ e 9' on ' n Potential derivative plane, £ = -(d<p/dX) + ifdyr/dY). 

f 0160 *' 31- Figure 8. - Auxiliary planes used in mapping procedure. 

Figure 7. - Porous cooled rectangular duct. 
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(b) Duct aspect ratio, 2. 

Figure 9. - Configurations of porous region for various duct aspect ratios and values of imposed physical 
parameter. 











Physical parameter, — L — ^ 

Figure 10. - Dimensionless heat and coolant mass flow rates at interior 
surface of porous duct as function of imposed physical parameter. 


NASA-Langley, 1973 




